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I. INTRODUCTION 

The existence of big bang singularity in the very early universe remains an open issue although 
most of the mysteries of standard hot big bang cosmological model have been addressed by inflationary 
epoch [1], To overcome this initial singularity, Ellis et al proposed a scenario called Emergent Universe, 
which is ever existing without singularity and has an almost static behaviour in the infinite past [2, 3]. 
Eventually the model evolves into an inflationary phase. 

Recently there have been lot of interest in emergent universe models based on standard as well as 
modified gravity [4-7] . This resurgence of interest is because of CMB observations favouring early 
inflationary universe [8, 9] and probably this is the possible mechanism for present acceleration [10]. 
Historically the credit goes to Harrison [11] for obtaining a model of closed universe with radiation 
and he showed that asymptotically it approached Einstein static universe (ESU). Then after a long 
gap Ellis and Maartens [2, 3] in recent past were able to formulate closed universe with a minimally 
coupled scalar field 4>- However, exact analytic solutions were not presented in their work and only 
asymptotic behaviour agreed with emergent universe. Subsequently, Mukherjee et al [12] obtained 
solution for Starobinsky model with features of an emergent universe. Also Mukherjee et al [13] 
formulated a general framework for an emergent universe using an adhoc equation of state which has 
exotic behavior in some cases. Very recently the idea of quantum tunnelling has been used to model 
emergent universe [14], 

In most of the General Relativity (GR) based models, Emergent universe is generally obtained 
for spatially closed past ESU which is not stable . Contrary to this in modified gravity, emergent 
scenario can be obtained in spatially flat as well as open/closed model of universe [4, 5]. Also in the 
context of modified gravity, including Braneworld models it has been found, a stable ESU can be 
obtained [15]. It may be noted that a stable ESU also exists in an open universe in f(T) gravity , loop 
quantum cosmology and Horova Lifshitz gravity [16-18]. Since astronomical observations including 
very recent Planck results favour spatially flat universe it is imperative to examine Emergent scenario 
in spatially flat universe. In this paper, we look for a possibility of emergent scenario from particle 
creation mechanism in spatially flat braneworld models. The modern era of brane cosmology began 
in the context of large extra dimensions, made possible by the hypothesis that the standard model 
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of particle physics is localized on a D-brane. Braneworld models of the universe have been the focus 
of attention in the current decade for the search of many outstanding problems in Cosmology. These 
kind of models are inspired by String Theory. In braneworld scenario, our four dimensional universe 
(a brane) is a hypersurface embedded in higher dimensional bulk spacetime. In this brane-bulk 
scenario, all matter and gauge interactions (described by open strings) are localised on a brane while 
gravity (described by closed strings) may propagate into whole space time. This means that gravity 
is fundamentally a higher dimensional interaction and we only see the effective 4D theory on brane. 
Among the different proposals for brane models, the two prominent are RS (proposed by Randall and 
Sundrum) and DGP(proposed by Dvali,Gabadaze and Porrati) models. In RS models, the hierarchy 
problem could be solved by a warped or curved extra dimension showing that fundamental scale could 
be brought down from the Planck scale to 100 GeV. 

While RS 2 model have only one (1 + 3)-brane [20], the RSI model have two two (1 + 3)-branes 
at the ends of the orbifold S 1 jZ<i [19]. Due to its simple and rich conceptual base, RS2 model is very 
popular and got much attention [21-24] and in particular, on the inflationary scenario. 

The main idea of the DGP model is the inclusion of a four dimensional Ricci-Scalar into the action 
.On the 4-dimensional brane the action of gravity is proportional to Mj, whereas in the bulk it is 
proportional to the corresponding quantity in 5-dimensions. The model is then characterized by a 
cross over length scale 

Ml 

Tc ~ 2M| 

such that gravity is 4-dimensional theory at scales a <C r c where matter behaves as pressure less 
dust but gravity leaks out into the bulk at scales a ^$> r c and matter approaches the behaviour of a 
cosmological constant [25-27]. 

In the standard cosmology RS2 model modifies the early universe where as DGP gravity modifies 
the late universe. So in the present universe brane corrections are not effective in RS2 model. But 
as phantom energy increases with the expansion, it brings drastic changes in the RS2 model [28]. 
This interesting and recent feature of RS2 model which can modify late time cosmic expansion (if the 
energy density of the matter content increases at late time e.g., phantom scalar field) has also been 
beautifully shown by using dynamical system tools in ref [29]. But this kind of modification does not 
appear if the energy density of the matter content decreases with cosmic expansion e.g, quintessence 
scalar held, radiation, dust etc. Thus RS brane model can also modify late time cosmic expansion in 
addition to its appreciable impact on early universe cosmology. 

Very recently Chakraborty in ref [6] formulated a emergent universe model in GR using the mech¬ 
anism of particle creation. The aim of this paper is to formulate an Emergent Scenario in Braneworld 
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as consequence a of particle mechanism. In a sense this paper is a generalisation of the one reported 
in [6] to include higher dimensional behaviour. Often using BG one obtains cosmological surprises. 
It may be noted that from thermodynamical aspect it has been proposed that entropy consideration 
favors the Einstein static phase as the initial state of our universe [30, 31]. The paper is organized 
as follows : Section 2 deals with non-equilibrium thermodynamics from the perspective of particle 
creation in cosmology while in section 3 we present Emergent scenario in Braneworld models. The 
paper ends with a short discussion and Concluding Remarks 4. 

II. NON EQUILIBRIUM THERMODYNAMICS AND PARTICLE CREATION IN 

COSMOLOGY 

Suppose E be the internal energy of a closed thermodynamical system having N particles. The 
first law of thermodynamics which is essentially the conservation of internal energy is given by [32] 

dE = dQ — pdV (1) 

where as usual p is thermodynamic pressure, V any co-moving volume and dQ represents the heat 
received by the system in time dt. 

From the above conservation equation, we get the Gibb’s equation given by 

Tds = dq = d(^—^J + p d(^—^J (2) 

E 

n = y , dq = jj- and p = — being the particle number density , heat per particle and energy 
density of the system respectively. It may be noted that this equation is also true when the particle 
number is not conserved i.e., the system is not a closed system [33]. Here we shall consider an open 
thermodynamical system where the number of fluid particles is not conserved and the non conservation 
of fluid particles is expressed as (Nfa / 0) 


h + Qn = nT (3) 

where = nu 11 is the particle flow vector, is the particle four velocity, 0 = is the fluid 
expansion, T stands for the rate of change of the number of particles (N = na 3 ) in a co-moving 
volume V = a 3 and by notation n = n 

The positivity of T indicates creation of particles while there is annihilation of particles for negative 
r . Any non zero T will behave as an effective bulk pressure of thermodynamical fluid and one can 
use non equilibrium thermodynamics [34]. 
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In this work, we shall consider a flat, homogeneous and isotropic model of the Universe, given by 
the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric in comoving coordinates (f, r, 9, < f>) as 

ds 2 = dt 2 — a 2 (t ) [dr 2 + r 2 {d0 2 + sin 2 9d(j) 2 )] , (4) 

where a(t) is the scale factor of the Universe. Here we consider universe as an open thermodynamical 
system which is non equilibrium in nature due to particle creation. 


The energy momentum tensor of the cosmic fluid is characterized by the energy momentum tensor 

T^u = (P + P + Ttyu^uv + {p + n )g fll/ (5) 

The energy conservation relation of T-y = 0 takes the following form 

p + 3H(p + p + n) = 0 (6) 

where H = — is the Hubble parameter. Here the pressure term n is related to some dissipative 
a 

phenomena (say bulk viscosity). In the present context, however the cosmic fluid may be considered 
as perfect fluid where dissipative term n is the effective bulk pressure due to particle creation or 
equivalently the conventional dissipative fluid is not taken as cosmic substratum , rather a perfect 
fluid with varying particle number is considered. For an isentropic (or adiabatic ) particle production 
this equivalence can be nicely described as follows [34, 35]: 

Using the conservation eqs. (3) and (6) the entropy variation can be obtained from Gibb’s equation 
(2) as 


nTs = — 3HU — T(p + p) (7) 

where T is the temperature of the fluid. 

But due to isentropic (or adiabatic ) process, the equilibrium entropy per particle does not change 
(as it does in dissipative process), i.e, s = 0, and from eq.(7) the effective bulk pressure is determined 
by particle creation rate as 

U = ~^H {P + P) (8) 

Hence, the bulk viscous pressure is entirely determined by particle production rate. So we may say 
that a dissipative fluid is equivalent to a perfect fluid having varying particle number. It may be 
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noted that, although s = 0, still there is entropy production due to enlargement of the phase space 
of the system. In the present context the same effect is due to the expansion of the universe. This 
effective bulk pressure does not correspond to conventional non equilibrium phase, rather a state 
having equilibrium properties as well. We note that it is not the equilibrium era with T = 0 

From eq.(8), we see that if we know the Hubble parameter H, then the particle creation rate can 
be determined. Further using Friedman equations one can relate T to the evolution of universe. In 
what follows as in [34] , assuming the particle creation rate as a function of the Hubble parameter , we 
can study the cosmological evolution in Braneworld. 


III. EMERGENT UNIVERSE IN BRANEWORLD SCENARIO 

A. RS2 model 

In a homogeneous and isotropic model of universe given by FRLW metric (4), the modified Fried¬ 
mann equation of RS2 model is given by 


H 2 



l + e 


2AJ 


(9) 


where e = ±1 and the sign of e are related to positive and negative brane tension respectively. (For 
simplicity we are using 8irG = 1 ) 

In the absence of bulk pressure, the matter content sector satisfy 


p + 3H(p + p) = 0, (10) 

In order to highlight some characteristics of this model, we review our paper [28] for the case e = — 1 
with p = puj , n = 0 and uj < —1. Using eqs. (9) and (10), we see that exact solution of energy density 
is given by 


P(t) 



( 11 ) 


where po < 2A is the current energy density of dark energy and to is the present time. From above 
equation we see that the braneworld model with e = 1 (A > 0) that describes an accelerated expansion 
ends in big-rip singularity in time 
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t 8 — ^0 


\/3|l + 


w\ 


1 


1 


2A + V Po + 2AJ 


( 12 ) 


as obtained in most of the GR based models. Brane gravity corrections make changes in the behaviour 
of phantom fluid in case phantom fluid dominates RS2 model of the universe with e = — 1 (A < 0). In 
this case it is found that it explains the present cosmic acceleration and exhibits deceleration after a 
finite time, when energy density grows sufficiently [28] . 

In the present work we take II / 0, the eq. (9) can be rewritten by introduction of dimensional 
variable x = p/ 2A 


3 H 2 = p eS (13) 

where the effective density is given by 

p e s = 2Ax(l + ex) (14) 

It is worth noting that when x « 1, the GR based cosmology is recovered and brane effects give new 
kind of behaviour when x » 1. Acceleration and deceleration of this model without bulk pressure 
have been studied in [28, 37] . The case of e = —1 is (also known as bouncing universe) have been 
extensively studied in [38, 39]. Also the effect of bulk brane viscosity in this model in the framework 
of Eckart’s theory has been studied in [40]. 

Taking derivative of eq. (13) w.r.t cosmic time t and using (6) , we get 


H = — - [p e ff + PeS + n eff ] 

where effective pressures are given by 

Peg = 2X[ux(l + e2x) + ex 2 ] 


(15) 


(16) 


n cff = n(i + e2x) 

and the effective equation of state is given by 


(17) 


Pei f T n eff 
PeS 


1 

1 + ex 


(1 + 2ex) 



+ ex 


k-’eff — 


(18) 
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FIG. 1. Variation of T with H in RS-II Model, considering e = 0.787 and / = —6.804. 


Adding eqs. (14) , (16) and (17) and using (8) , we can rewrite eq (15) as 


2 H = -(l + e2®)(l+w)2Ax 1 


m) 


Again rewriting modified Friedmann equation (13) in terms of x, we get 


2e\x = e\[B RS (e,H) - 1] where B RS (e,H) = 


6eH 2 


+ 1 


Substituting the value of 2Ax in eq. (19), we get 


2 H = 


6(1 + uj)H 2 


[1 + 1/ B R s(e, H)] V 3 H) 

As the particle creation rate is a function of Hubble parameter H [34] , we choose 


(19) 


( 20 ) 


( 21 ) 


r = 3H 


f 1 f 

(1 — e) + —-(e — —) 

V ^ H B RS y H ’J 


( 22 ) 


where e and / are constants. 

Thus using eqs. (21) and (22), we get the differential equation of H describing cosmic evolution as 


2H = —6(1 + u)(eH 2 — fH) 


(23) 


So the Hubble parameter can be obtained by solving eqs. (23), as (w ^ — 1), 

Ho 


H 


= H i + exp 




(24) 
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where H 0 = 6(1 + u)f , Hi = 6(1 + u)e and t 0 is the constant of integration. Hence the scale 
factor can be obtained from the above expression for Hubble parameter as 


El 

2 


= 1 + Hi exp 


H 0 (t - t 0 ) 


(25) 


where a o = a(to) 

The above cosmological solution shows the following asymptotic behavior for / < 0: 

(i) a —>• a-o , H —>• 0 as t —>• —oo 

(ii) a ~ ao , H ~ 0 for t « to 

{in) a ~ exp[fft(i - t 0 )} , H ~ for t » t 0 • 

The above asymptotic features show that the above cosmological solution (described by equations 
(24) and (25) ) describes a scenario of emergent universe. 


B. DGP model 

The modified Friedmann equation in this case is given by 

H 2 -e- = P - (26) 

r c 3 

M 2 

where r c = ttttz is the crossover scale which determines the transition from 4D to 5D behavior and 

c 2 Mg 

e = ±1. For e = 1, we have standard DGP(+) model which is self accelerating model without any form 
of dark energy, and effective oj is always non phantom. However for e = —1, we have DGP(-) model 
which does not self accelerate but requires dark energy on the brane. It experiences 5D gravitational 
modifications to its dynamics which effectively screen dark energy. 

The above equation can be rewritten as 

3H 

3 H 2 = p e ff where effective density is p e g- = p + e- (27) 

r c 

Taking derivative of eq. (26) w.r.t cosmic time t and using eq. (6) , we get 


(2 H - e/r c )H = -H(p + p + U) 


Now using eq. (8) in above equation we get 

2 H = —6H 2 (1 + w)(l - —^jBv GP {e,H) 


where 


-Bdgp(g H ) 


H - e/r c 
2 H - e/r c 


and it gives brane effects 


(28) 


(29) 


(30) 
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_ H _ 

- c = -2.5 , d = 6.673 - c = -4.5 , d = 8.673 

■ ■ ■ ■ c = -6.5 , d = 10.673-c = -8.5 , d = 12.67: 


FIG. 2. Variation of T with H in DGP Model, considering e = —1 and r c = 4 . 

In this case , we choose the particle creation rate as 

T = 3H[1 — (c —d/H)/BDQp(e, H)], c and d are constants (31) 

As before using modified Friedmann eqs. (29) and (31) we get the differential equation of H 
describing cosmic evolution of DGP model as 


2H = -6(1 +u)cH 2 + 6d(l +u)H (32) 

This evolution equation is identical in form to that for RS-II brane gravity (i.e. eq.(23) ) and hence 
it is also poosible to have emergent scenario in DGP brane model in the framework of particle creation 
mechanism. 


IV. SHORT DISCUSSION AND CONCLUDING REMARKS 

The present work deals with spatially flat braneworld models in non-equilibrium thermodynamics, 
based on particle creation mechanism. For simplicity, we assume the thermal process to be adiabatic 
in nature and as a result the particle creation rate is related to the cosmic evolution. In both RS-II and 
DGP brane model, cosmic fluid with constant equation of state (p = up , w / — 1) is chosen. In both 
the brane models the particle creation rate is chosen as different functions of the Hubble parameter 
and the brane effect is incorporated through the functions Bp s and Rdgp respectively. However the 
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cosmic evolution equation is similar in both the models and we have emergent scenario in the early 
phase of the evolution. Further it should be noted that the present model of emergent scenario due 
to particle creation is on ever expanding model of the universe. As the particle creation rate is not 
constant, so the present model has the basic difference with the steady state theory of Fred Hoyle et 
al [41, 42] in contrast to the emergent scenario in Einstein gravity with constant particle creation rate 

[6]. Moreover, one can easily check that if we remove the brane effect (i.e. A —>• oo in RS2 or r c —>• oo 
in DGP model) then the particle creation rate may be taken to be constant (by choosing e = 1/2 
in RS2 model and c = 1/2 in DGP brane model). The variation of T with respect to the Hubble 
parameter has been presented in figs. (1) and (2) for RS2 and DGP brane model respectively for 
various choices of the parameters involved. Finally, we conclude that it is possible to have emergence 
scenario in brane gravity models for non-equilibrium thermodynamics in the framework of particle 
creation mechanism. 
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